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Energy levels of the four lowest bands in I60,i62,i64j3y ^^^^ lesgj.^ B{E2) 
transition strengths between the levels, and the B(M1) strength distribution 
of the ground state, all calculated within the framework of pseudo-SJ/a 
model, are presented [9] . Realistic single-particle energies and quadrupole- 
quadrupole and pairing interaction strengths fixed from systematics were 
used in the calculations [6]. The strengths of four rotor-like terms, all 
small relative to the other terms in the interaction, were adjusted to give 
an overall best fit to the energy spectra. The procedure yielded consistent 
parameter sets for the four nuclei. 
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1. Introduction 

Experimental nuclear physicists, with improved equipment and tech- 
niques, continue to challenge nuclear theorists with interesting new phe- 
nomena. The measurement of new levels, some lying below 3 MeV, raises 
questions about the nature of collective excitations in atomic nuclei. Even 
though heavy deformed nuclei with ^ > 150 are good candidates for prob- 
ing such degrees of freedom, microscopic calculations are very important 
to gaining a deeper understanding of the nature, for example, of low-lying 
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excited O"*" bands and K'^ = 4+ states. Of special interest to us is the 
fragmentation of the ground state Ml strength distribution. 

Over the last several years, various properties of low-energy states in 
^6°Dy [1], i^^Dy [2], ^^^Dy [3] and ^^SEr [3, 4, 5] as well as other rare earth 
nuclei have been measured. In a recent article, we presented a description of 
the even-even 156,158,160q^j isotopes [6] within the framework of the pseudo- 
SU3 model. One of our main goals of the present study was to continue to 
test the Hamiltonian and obtain a consistent set of parameters for a larger 
set of nuclei. In order to achieve this goal we fixed the strengths of the 
dominant terms in the interaction: single-particle energies as well as the 
two-body quadrupolc-quadrupole and pairing interactions. Four smaller 
terms were varied to "fine tune" the results to experiment. 

In this paper we apply the pseudo-St/s model to the normal parity bands 
in the I60'i62,i64j-)y ^^^'En, These nuclei, as for the Gd isotopes studied 
earlier, all exhibit well-developed rotational ground state bands as well as 
states that are associated with low-lying K'^ = 0+ and K'^ = 2+ bands. 
Here wc provide a theoretical description for states of the low-lying bands 
in each of these isotopes including B(E2) transition strengths between the 
states and the ground state B(M1) strength distribution. As appropriate, 
we also give results for the = 4+ and = 1+ bands in some of 
these nuclei. The theory yields a reasonable reproduction of the observed 
low-lying 1+ states, the ground state Ml sumrule and its energy-weighted 
centroid, as well as the observed fragmentation of the Ml strength but not 
to the same level of accuracy as for other properties. 

The Hamiltonian we used in the study is given in Section 2. Excitation 
energies, intra-band B(E2) transition strengths, and ground state B(M1) 
strength distributions are reported in Section 3. Some conclusions that 
follow from the analyses are offered in Section 4. 

2. Model space and parameters 

The nuclei considered in the present calculations have closed shells at 
JYtt = 50 for protons and N^, = 82 for neutrons. To build basis states, we 
considered an open oscillator shell for each (% = 4 for protons and rj^ = 5 
for neutrons) along with their intruder state complements (/111/2 for pro- 
tons and ii3/2 for neutrons) even though particles in these unique-parity 
intruder levels were only considered to renormalize the normal-parity con- 
figurations through the use of an effective charge. These oscillator shells 
have a complementary pseudo-harmonic oscillator shell structure given by 
ifcr {a = TTjiy) = Tjcr — 1- Approximately 20 pseudo-Sf/a irreducible repre- 
sentation (irreps) with largest values of the second order Casimir operator 
C2 {Q ■ Q = 4C2 — 3L^) were used in building the basis states. 



draayer printed on February 8, 2008 



3 



We used the following realistic pseudo-5J73 Hamiltonian: 

H = + Hgp - - X Q ■ Q - G^TT Hp - Gv Hp 

+ a + h K'^j + a^Cz + asym C2 . (1) 

Strengths of the quadrupole-quadrupole {Q ■ Q) and pairing interactions 
(Hp) were fixed, respectively, at values typical of those used by other au- 
thors, namely, x = 35 A^/^ MeV, = 21/ A MeV and G^, = 19/ A MeV. 
The spherical single-particle terms in this expression have the form 

H^p = J2(CaK-^i. + D^t)- (2) 

Since only pseudo-spin zero states were considered, matrix elements of the 
spin-orbit part of this interaction vanish identically. Calculations were car- 
ried out under the assumption that the single-particle orbit-orbit (P) inter- 
action strengths were fixed by systematics [7] , 

D„{a = TT, Z/) = TlUKajJia , (3) 

where Tiw = 41/ A^/^ with and assigned their usual oscillator values [7]: 

K„ = 0.0637, /i^ = 0.60 

= 0.0637, = 0.42 . (4) 



parameter 
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7.55 


X X 10-3 


6.84 


7.12 


7.27 


7.42 




-0.283 


-0.286 


-0.287 


-0.289 




-0.198 


-0.200 


-0.201 


-0.202 




0.125 


0.128 


0.130 


0.131 


Gy 


0.101 


0.104 


0.105 


0.106 


a xlO-3 


-2.1 


-2.0 


0.0 


1.0 


b 


0.022 


0.00 


0.08 


0.10 


O'sym X 10 


0.80 


1.20 


1.40 


1.45 


03 X lO"'' 


0.75 


0.65 


1.32 


1.36 



Table 1. Parameters of the pseudo-SiJa Hamiltonian. 



Relative excitation energies for states with angular momentum 0"*" are 
determined mainly by the quadrupole-quadrupole interaction. The single- 
particle terms and pairing interactions mix these states. With the strength 
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of these interactions fixed as in Table 1, the states he very close to their 
experimental counterparts while the O3" states usually slightly above the 
experimental ones. Of the four "free" parameters in the Hamiltonian, a was 
adjusted to reproduce the moment of inertia of the ground state band, 03 was 
varied to yield a best fit to the energy of the second 0"*" state (the energy of 
the third 0+ was not included in the fitting and as the results given below 
show these all fall slightly higher than their experimental counterparts), 
O'sym was adjusted to give a best fit to the first 1+ state, and b was fit to 
the value of the band-head energy of the = 2+ band. 

In the rotational model the projection K of angular momentum on 
the body-fixed symmetry axis is a good quantum number. For each in- 
trinsic state with a given value of K there is a set of levels with L = 
K, K + 1, K + 2, . . . ., except for K = when L is either even or odd 
depending on the intrinsic {D2) symmetry of the configuration. Elliott 
[8] used group-theoretical methods to investigate classification schemes for 
particles in a three-dimensional harmonic oscillator potential for which the 
underlying symmetry is 5*^73. He noted that the angular momenta in an 
irrep of 5*^73 can be grouped in a similar way to that of the rotor, the differ- 
ences being that there are a fixed number of K values and that each band 
supports a finite number of L values rather than being of infinite length. 
The angular momentum content of an SU3 irrep (A, /i) can be sorted into 
K bands according to the following rule [9]: 

K = imii{X, fj,) , min(A,/x)— 2, . . . , 1 or 0, (5) 

where 

L = (A + m), (A + m) - 2, . . . , 1 or (6) 

for iiT = and 

L = K, K + 1, K + 2, ...,{X + ,j,)-K + l (7) 

for K ^ 0. Hence, the allowed angular momentum values in the leading 
SU3 irrep for ^^"Dy with (28,8) are L = 0, 2, ...,36 for the if = band, 
L= 2, 3, 35 for the K = 2 band, etc. 

3. Results 

Using the Hamiltonian and procedure described in the previous section, 
good agreement was obtained between the experimental and calculated en- 
ergies of the first three low-energy bands in each of the nuclei considered. 
States in the ground-state, first excited K = 2 and first excited K = {) bands 
were all found to lie very close to their experimental counterparts. A second 
excited K'^ = 0"*" band was identified at approximately 0.5 MeV above its 
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Energy spectra of ^^°Dy, ^''^Dy, and ^^^Dy obtained using Hamiltonian 1. 
represents the experimental results and 'Th.' the calculated ones. Figures 
and f) give the theoretical and experimental Ml transition strengths from 
= ground state to the various J = 1 states. 



experimental counterpart. As noted above, this level was not included in 
the fitting procedure. Two other bands, = 1'^ and K'^ = 4+, were also 
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identified for eacli of the nuclei. 
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Fig. 2. Energy spectra of ^^^Er obtained using Hamiltonian (1). 'Exp.' represents 

the experimental results and 'Th.' the calculated ones. Figures b) gives the theo- 
retical predictions for Ml transition strengths from the J = ground state to the 
various J = 1 states. 



As can be seen from Table 1, the parameters that the fitting procedure 
yielded for the Dy isotopes vary smoothly from one nucleus to another. In 
fact all four of the 'free' parameters decrease as the mass number A increases. 
Results were also generated for ^^^Er which has one pair of protons more 
than the Dy isotopes. Strengths of the interactions in this case are given 
in second column of Table 1. Since A is larger for ^^^Er than for the Dy 
isotopes, one might expected even smaller values for the parameters, and 
while this is what is observed for a and asym, the extra proton pair requires 
values for b and as that are smaller than the ones for I60,i62£)y larger 
than those for ^^^Dy. 

Figure 1(a) shows the calculated and experimental [10] K = 0, K = 2 as 
well as the first and second excited K = bands for ^^"^Dy. The agreement 
between theory and experiment is excellent for the first three bands with 
relative differences between calculated and experimental energies being in 
all cases less than 7%. The model predicts a continuation of the first excited 
K = band with two additional states of angular momentum 6 and 8 in 
^^'^Dy. The second excited K = band is identified to be about 0.5 MeV 
higher than the experimental one. As for ^^'^Dy, the agreement between 
theory and experiment is excellent for the first three bands in ^^^Dy and 
^^^Dy. Again, the second excited K = band is identified at slightly higher 
than the experimental one. The experimental and calculated energy spectra 
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for these nuclei are given in Figs. 1(c) and 1(e). 

The energy spectra for ^^^Er is given in Fig. 2 (a). Even though this 
nucleus has one additional pair of protons, the calculated energy spectra 
exhibits the same behavior as is Dy isotopes. The ground-state, K = 2, and 
first excited K = bands are well reproduced. An angular momentum 8 
state is identified as a continuation of the first excited K = band. The sec- 
ond excited K = band is identified and lies higher than the experimental 
one. As Fig. 1 shows, the difference in energy between the calculated and 
experimental band heads of the second excited K = bands increases as 
the mass number A increases. For ^^''Dy and ^®^Dy the difference is about 
0.3 MeV while for ^^^Dy it is about 0.5 MeV, and for ^^SEr is about 0.65 
MeV. The states in the second excited K = band are well defined. The 
results suggest that it may be interesting to see how the parameters change 
with changing proton rather than neutron number. 

Collective models have emerged as a result of attempts to reproduce 
experimental observations with simple calculations. In the evolution of the 
SUs model, and later the psendo-SUs model, one of the main motivations 
was to achieve a good description of deformed nuclei using a small configu- 
ration space. In the present study the configuration spaces are small com- 
pared with typical ones that are used in shell-model calculations based on 
m-scheme configurations. To provide a better understanding of the theory 
it is useful to take a closer look at the eigenvectors of the states in the four 
bands under consideration. This is done in Table 2 where the SU3 content 
of the calculated eigenvectors for states of the four lowest bands in ^^^Dy 
are given. The percentage distributions of each eigenvector across the (A, /x) 
values are given in the second column. Only basis states that contribute 
more than 2% are identified. In the ground state band there is clearly a 
dominant SU^ irrep which is the leading irrep in the coupling of the leading 
proton and neutron SU3 irreps. An additional five irreps contribute in total 
with about 40 % to the eigenvector strength. The same leading irrep is the 
most important in the K = 2 band, and as we will see later in the K = 4 
band. Again there are only about six irreps that contribute more than about 
2% to the eigenvectors. The first excited K = band exhibits a different 
dominant SU3 irrep, namely 60% [(4, 10)^(18, 4)^] (22, 14) for ^^^Dj. The 
second excited K = state is a strong mixture of two S'f/s-irreps, 47% of 
[(10, 4)^ (20,0)^] (30, 4) and 34% of [(10, 4)^(18, 4)^] (28, 8). As before, only 
a couple irreps contribute to more than about 2% to the eigenvectors. 

The SU:i content of the eigenvectors is fairly constant across the states 
within a band. The percentages vary slowly and smoothly from that of the 
band-head as one moves up the band to states with higher values of the 
angular momentum. To see this, the SU3 content of calculated eigenvectors 
for the states in the ground-state band in ^^^Dy are given in Table 3. Only 
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Table 2. SUz content of calculated eigenvectors for states of the four low- lying 
band-heads in ^^^Dy. The percentage distributions of each eigenvector across the 
(A, /x) values are given in the second column. Only basis states that contribute to 
more than 2% are identified. 
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Table 3. SUz content of calculated eigenvectors for the states in the ground-state 
band in ^^^Dy using parameters from Table 1. Only the basis states that contribute 
to more than 2% are identified. 
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3.1. B(E2) transitions 

Theoretical and experimental [10] B(E2) transitions strengths between 
the states in the ground state band in -"^^^Dy are shown in Table 4. The 
agreement between the calculated and experimental numbers is excellent. 
The B{E2;2i 4i) is equal to within 1% of the experimental value, and 
the last two calculated B{E2) values differ from the experimental values by 
less than 0.1 e^6^ which is well within the experimental error. Excellent 
agreement with experimental B{E2) data is also observed in ^^^Dy and 
I64j])y Contributions to the quadrupole moments from the nucleons in the 
unique parity orbitals are parameterized through an effective charge [9], ey, 
with Ci, = Sf, and = 1 + e/, so the E2 operator is given by [9] 

Theoretical intra-band B(E2) transition strengths between the states in the 
if = 2 as well as the first and second excited K = Q bands are given 
in Table 5. Note that the strengths of the transitions probabilities are 
consistent across all four bands (Tables 4 and 5). 



Ji 




B{E2- Ji ^ Jf 


) {e^h') 






Exp. 


Theory 


Oi 


^2i 


5.134 ± 0.155 


5.134 


2i 


^4i 


2.675 ± 0.102 


2.635 


4i 




2.236 ± 0.127 


2.325 


6i 


^8i 


2.341 ± 0.115 


2.201 



Table 4. Experimental and theoretical B{E2) transition strengths between mem- 
bers of ground state band of ^^^Dy. 



3.2. Ml transitions and the K = 1]^ hand 

Another test for the theory is the Ml transition strength distributions 
that can be obtained using eigenvectors of the diagonalized Hamiltonian 
(1). The calculated and experimental Ml strength distributions for the Dy 
nuclei are given in Figs. 1 b), d) and f). The calculated Ml strength 
distribution for ^^^Er is given in Fig. 2 b). For illustrative purposes, the 
energies and Ml transition spectra are given opposite one another. 

The starting point for a geometric interpretation of the scissors mode 
within the framework of the SU^ shell model is the well-known relation of the 
SUz symmetry group to Rot{3), the symmetry group of the rotor [11, 12]. 
The structure of the intrinsic Hamiltonian allows for a rotor-model inter- 
pretation of the coupled SU^ irreps (X^jfJ-Tv) and (Ai^,/iiy) for protons and 
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Table 5. Theoretical B{E2) transition strengths between states of the K = 2, 

= O2, and K = O3 bands of '^^^Dy. The energies are labeled with the subindex 
7 for the K = 2 band, a, and b for the first and second excited K = Q bands. 



Nucleus 


Ei?(Ml)[/x^,] 


Experiment 


Calculated 


Pure SU^ 


Theory 




2.48 


4.24 


2.32 




3.29 


4.24 


2.29 


164Dy 


5.63 


4.36 


3.05 



Table 6. Total B(M1) strength from experiment [10] and the present calculation. 



neutrons, respectively. According to the Littlewood rules [13] for coupling 

Young diagrams, the allowed product configuration can be expressed in 
mathematical terms by using three integers (m, 1, k): 

(Att, A^Tr) ® (Ai/, iJiu) = ®m,l,k{^7T + \u - "^rn + I, (XT, + fx„ - 21 + m)^, (9) 

where the parameters / and m are defined in a fixed range given by the 
values of the initial SU^ representations. In this formulation, k serves to 
distinguish between multiple occurrences of equivalent (A,//) irreps in the 
tensor product. The number of k values is equal to the outer multiplicity, 
Pmax {k = 1,2, .. . ,Pmax)- The I and m labels can be identified [15] with 
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excitation quanta of a two-dimensional oscillator involving relative rotations 
the angle between the principal axes of the proton and neutron system, 
and 0, the angle between semi-axes of the proton and neutron system) of 
the proton-neutron system, 

m = riQ, I = n^. (10) 

These correspond to two distinct type of l"*" motion, the scissors and twist 
modes, and their realization in terms of the pseudo-S'C/a model. 

The SU3 irreps obtained from the tensor product (9) that contain a. = 
1"*" state are those corresponding to {m,l,k) = (1,0,1), (0,1,1), (1,1,1), 
and (1, 1,2). A pure 5*^73 picture gives rise of a maximum of four l"*" states 
that are associated with the scissors, twist, and double degenerate scissors- 
plus-twist modes [(1,1,1) and (1,1,2)] [15]. Results for the Dy isotopes, 
assuming a pure pseudo-S'C/3 scheme, are given in Table 7. 

The experimental results [10] given in Figs. 1 b). d), and f) suggest a 
much larger number of l"*" states with non-zero Ml transition probabilities 
from the 0"*" ground state. The SU3 breaking residual interactions lead to 
a fragmentation in the Ml strength distribution, since the ground state 
0+ is in that case a combination of several SU3 irreps, each with allow 
Ml transitions to other SUs irreps. Overall, the total Ml strength is in 
reasonable agreement with the experimental results (Table 6). In ^^^Dy the 
total Ml strength is slightly underestimated, which may be due to spin 
admixtures in the wavefunction, which is not included in this work. 



Nucleus 


(Att, /^tt) 


(Att, /^tt) 




(A,/x)i+ 


B(M1) 


mode 


160-162-Qy 


(10,4) 


(18,4) 


(28,8) 


( 29, 6) 
( 26, 9) 
(27,7)1 
(27, 7)2 


0.56 
1.77 
1.82 
0.083 


t 
s 

s-l-t 
t-l-s 


iti^Dy 


(10,4) 


(20,4) 


(30,8) 


(31,6) 
(28,9) 
(29,7) 
(29,7) 


0.56 
1.83 
1.88 
0.09 


t 

s 

s+t 
t+s 



Table 7. B(M1) transition strengths in the pure symmetry limit of the pseudo 
SU^ model. The strong coupled pseudo-Sf/a irrcp (A, /i)^^ for the ground state is 
given with its proton and neutron sub- irreps and the irreps associated with the 1+ 
states, (A',/i')i+. In addition, each transition is labeled as a scissors (s) or twist 
(t) or combination mode. 



According to Eq. (5), and depending upon the values of the (A, /x) irreps, 
there are several K = 1 bands. Since we are interested in Ml transitions 
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Ji Jf 


B{E2- Ji ^ Jf) 


Energy(J/) [MeV] 
Exp. Th. 








1.746 


1.749 


ll 


^25 


1.312 


1.783 


1.797 


25 


^33 


0.944 


1.840 


1.905 


25 


^45 


1.170 


1.904 


1.941 


45 


^53 


0.802 




1.993 


45 


^64 


0.648 




2.066 


64 




1.291 




2.150 


64 


^84 


1.203 




2.292 



Tabic 8. The experimental and calculated energies of the K = \ band in -'^^^Dy 
using parameters from table 1. The B(E2) transition probabilities are given in the 
second column. The energy of the J = 1 state is given in the first row. 



from the ground state to the 1"*" states, several (A, ^) that allow K = 1 states 
are included in the configuration space. As a result, the model predicts 
several K = 1 bands. The first calculated K'^ = 1''" state lies very close to 
the first experimental one (the numbers arc identical in the first three digits) 
in ^^°Dy and ^^^Dy. Moreover this being the band-head of a = 1"*" band, 
the band is also remarkable well described. The calculated and experimental 
energies of the K'^ = l"*" band in ^^^Dy are given in Table 8 together with the 
B(E2) transition strengths. A few additional states of angular momentum 
5, 6, 7, and 8 are predicted as members of the first K = 1 band in ^^^Dy. 

3.3. K"" = 4+ band 

In the pure SU^ symmetry limit there are no interaction terms that mix 
states of different (A, /v,) irrcps. The energies of the nonzero band-heads 
associated with nonzero K values are then fixed, by-in-largc by the 6 pa- 
rameter (though band mixing within an irrep is allowed this is normally 
small, vanishing in the limit of large prolate configurations). A simple exer- 
cise yields a value of 2^ x 6 for the K = 2, J = 2 state and a value of 4^ x 6 
for the K = 4, J = 4 state. If, for example we assume the h value from Ta- 
ble 1 for ^^°Dy, the values for these two states would be 0.4 MeV for Jk=22 
and 1.6 MeV for Jk=44, respectively. However in the present version of the 
model these states are mixed due primarily to the single-particle energies 
and pairing interactions. In the present calculations b parameter adjusted 
to give a best fit to the band-head of the K = 2 state only, yielding a value 
of 2.29 MeV for the band-head of the K = A band in ^^^Dy, which is only 
slightly higher than the experimental value. 

In these nuclei there is an experimentally known low-lying = 4"*" 
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Ji Jf 


Energy(J/) [MeV] 
Exp. Th. 


B{E2; Ji ^ Jf) 
(e262) 




46 


1.536 


1.977 




46 


^54 


1.634 


2.168 


0.360 


46 


^66 


1.752 


2.209 


0.976 


66 


^74 


1.888 


2.514 


0.125 


66 


^86 




2.527 


1.784 



Table 9. Experimental and calculated energies of states of the angular momentum 
indicated by the right-side entry of the first column for members of the K = A 
band of ^^^Dy. The B(E2) transition probabilities are given in the last column. 
The energy of the J = 4 band-head is given in the first row. 



Ji 




Energy(J/) [MeV] 


B{E2-Ji^Jf) 




45 


2.173 1.836 




45 


^54 


1.944 


1.147 


45 


^67 


2.124 


0.493 



Table 10. Experimental and calculated energies of states of the indicated angular 
momentum given by the right-side entry of the first column for members of the 
K = 4 band in ^^^Dy. The B(E2) transition probabilities are given in the last 
column. The energy of the J = 4 band-head is given in the first row. 



band. A comparison between the experimental and calculated energies are 
given in Table 9 together with B(E2) transition probabilities for -"^^^Dy. One 
additional state of total angular momentum 8 is identified as a member of 
the K = 4: band, by both a strong B{E2) transition to the state of angular 
momentum 6 and through similar SU3 content with the other states in this 
band. In the ^^^Dy case there is an experimentally known K = 4 state at 
2.173MeV . The theory predicts a K = A state at lower energy (1.836Mey) 
and moreover according to the theory it should be a band-head. 

4. Conclusions 

This study of I60,i62,i64-Qy ^^Sgj, shows that pseudo-spin zero neu- 
tron and proton configuration with a relatively few pseudo-S'C/3 irreps with 
largest C2 values suffices to obtain good agreement with known experimental 
results. The Hamiltonian that was used included single-particle energies, the 
quadrupole-quadrupole interaction, and neutron and proton pairing terms, 
all with strengths fixed by systematics, plus four smaller rotor-like terms 
with strengths that were varied to maximize agreement with observations. 
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A consistent set of 'free' parameters was obtained. The results generated 
extended beyond quantities that were used in the fitting procedure, includ- 
ing intra-band B(E2) strengths, the Ml strength distribution of the ground 
state, and band-head energies of the first K'^ = 1^ and K'^ = 4+ bands. 

The Ml strength distributions were not fit to the data. Nevertheless, 
in all cases the summed strength was found to be in good agreement with 
the experiment numbers. The pseudo-S'f/3 model therefore offers a micro- 
scopic shell-model interpretation of the "scissors" mode [14], and in ad- 
dition, it reveals a "twist" degree of freedom that corresponds to allowed 
relative angular motion of the proton and/or neutron sub-distribution [15]. 
By adding one-body and two-body pairing interactions to the Hamiltonian, 
it was possible to describe the experimentally observed fragmentation of the 
Ml strength. The results suggest that more detailed microscopic descrip- 
tion of other properties of heavy deformed nuclei, such as g-factors and beta 
decay, may finally be within reach of a bona fide microscopic theory. 
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